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Abstract 

The optimization of antibodies is a desirable goal towards the development of better therapeutic strategies. The antibody 
1 1K2 was previously developed as a therapeutic tool for inflammatory diseases, and displays very high affinity (4.6 pM) for 
its antigen the chemokine MCP-1 (monocyte chemo-attractant protein-1). We have employed a virtual library of mutations 
of 11K2 to identify antibody variants of potentially higher affinity, and to establish benchmarks in the engineering of a 
mature therapeutic antibody. The most promising candidates identified in the virtual screening were examined by surface 
plasmon resonance to validate the computational predictions, and to characterize their binding affinity and key 
thermodynamic properties in detail. Only mutations in the light-chain of the antibody are effective at enhancing its affinity 
for the antigen in vitro, suggesting that the interaction surface of the heavy-chain (dominated by the hot-spot residue 
PhelOl) is not amenable to optimization. The single-mutation with the highest affinity is L-N31R (4.6-fold higher affinity 
than wild-type antibody). Importantly, all the single-mutations showing increase affinity incorporate a charged residue (Arg, 
Asp, or Glu). The characterization of the relevant thermodynamic parameters clarifies the energetic mechanism. Essentially, 
the formation of new electrostatic interactions early in the binding reaction coordinate (transition state or earlier) benefits 
the durability of the antibody-antigen complex. The combination of in silico calculations and thermodynamic analysis is an 
effective strategy to improve the affinity of a matured therapeutic antibody. 
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Introduction 

The exquisite specificity and high affinity of antibodies are 
increasingly exploited for therapeutic and biotechnological pur- 
poses, such as in cancer immunotherapy, in diagnosis, and in 
molecular bio-sensors [1—3]. Biobetter, defined as the biomedicine 
successfully developed through functional and/or physicochemical 
improvement of a natural molecule, provides a promising strategy 
for the next generation of therapeutics [4]. However, the 
optimization of the physicochemical properties of an antibody 
without perturbing the affinity and specificity for the cognate 
molecule is a challenging endeavor. Progress in this field is 
accelerating, and successful examples of improved antibodies using 
biobetter strategies have been reported, such as higher affinity for 
the cognate, increased stability in solution, enhanced pharmaco- 
kinetics, diminished immunogenicity, or conjugation to drug 
delivery systems [5-9]. Nevertheless, to advance further the 



deployment of biobetter strategies in the design and preparation of 
improved therapeutics it is necessary to strengthen our under- 
standing of the physicochemical properties of engineered antibod- 
ies [10]. 

The interaction surface of antibody-antigen complexes involves 
a large number of residues and water molecules establishing 
multiple non-covalent interactions that are difficult to quantify at 
the molecular and atomic levels. Although the evaluation of 
protein-protein binding energy remains a challenging task, 
progress in computing performance and force-field parametriza- 
tion are rapidly advancing our predictive capabilities [11,12]. 
Thus structure-based computational techniques are increasingly 
employed in the design of biotherapeutic antibodies [13]. 

Herein, we report the improvement of the binding affinity of the 
mature antibody 11K2 for its target antigen MCP-1 (monocyte 
chemotactic protein-1), an important therapeutic target in 
inflammatory diseases such as arteriosclerosis [14], allergy [15], 
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and rheumatoid arthritis [16]. The design of an optimized version 
of the antibody 11K2 with enhanced binding capabilities may 
improve its therapeutic value in the treatment of inflammatory 
diseases. Although the binding affinity of 11K2 for its antigen is 
very high (4.6 pM) [17], the analysis of the crystal structure of the 
antibody-antigen complex (PDB entry code 2BDN) suggests how 
to further increase the affinity by optimizing the complementary 
determining region (CDR) of the variable domain of the light 
chain (Vl) (Figure 1) [18]. Indeed, the interaction surface between 
the CDR of the V L of 1 1K2 and the antigen is 257 A 2 , a value 
significandy smaller than that of the CDR of the variable domain 
of the heavy chain (Vh) (51 1 A 2 ) [19]. Moreover, the evaluation of 
the thermodynamic properties of the engineered antibodies is a 
desirable goal towards increasing our predictive capabilities in 
future optimization strategies. 

In this report we have calculated the binding energy between 
MCP-1 and a virtual library of systematic mutations at each of the 
62 residues of the CDR of the antibody 11K2 by computational 
methods. The in silico calculations suggest that single-mutations 
carrying charged residues at specific locations favor the binding 
energy with the antigen. The affinity of the most promising 
candidates was determined experimentally by the technique of 
surface plasmon resonance (SPR), which identified several muteins 
with enhanced affinity for the antigen. The thermodynamic 
analysis reveals the fundamental mechanism explaining the 
superior binding capabilities of the optimized antibody. Our study 
highlights the benefits of combining in silico and in vitro 
methodologies for a more effective molecular design by biobetter 
strategies. 

Materials and Methods 

Comprehensive mutagenesis in silico 

Each of the 62 residues belonging to the CDR of 11K2 (as 
defined by Kabat [20] and Chothia [21]) were subjected to 
systematic mutagenesis in silico with each of the other 19 natural 
amino acids (62x19 = 1,178 mutations). The initial coordinates 
of 11K2 in complex with MCP-1 were retrieved from the PDB 
(entry code 2BDN). For each mutation in the antibody, 100 
randomized models were generated using the default parameters 
of the program MODELLER of the Discovery Studio Suite 
(Accerlys, San Diego, CA) [22]. Each model of the mutated 
antibody-antigen complex was optimized by a combination of 
simulated annealing followed by molecular mechanics minimiza- 
tion until the mean square gradient decreased below 0.01 kcal 
mol 1 A 2 . The interaction energy between antigen and antibody 
in each of the 1 1 7,800 models was calculated using the default 
values in the module Einteract of the package software MOE 
(Chemical Computing Group, Canada) using the AMBER99 force 
field after hydrogen atoms were explicitly added and minimized. 
We used the default dielectric constant (80.0), mimicking the 
behavior of water. Structural models of muteins lacking non- 
covalent interactions with the antigen were not considered for 
further analysis. For each mutation tested, the electrostatic and 
van der Waals interaction energies were summed up and applied 
as a surrogate to evaluate the affinity of the mutant for the antigen, 
as described previously [23]. As a reference we built 1,000 model 
structures of native antibody by optimizing the x-ray structure 
using the program MODELLER as explained above. Mutations 
displaying histograms with more favorable energy of interaction 
than that of wild type protein were selected for in vitro 
examination (H-L27R, H-L27K, H-N28D, H-N28Q, H-D31E, 
L-Y30K, L-N31R, L-N31K, L-S53D, L-S53E, L-T56D, L-T56E) 
(see a representative example in Figure SI). 



Expression and purification of MCP-1 

The DNA encoding MCP-1 was synthesized by GeneArt 
(Regensbyrg, Germany) and sequence-optimized for expression in 
Escherichia coli (Table SI). The MCP-1 gene was expressed in a 
vector pET26b (Novagen) displaying a hexa-histidine tag at the C- 
terminus. The DNA sequence was confirmed by the dideoxy 
chain-termination method. 

E. coli strain Rosetta2 (DE3) cells (Novagen) transformed with 
the expression vector of MCP-1 were grown at 28°C in 2x YT 
broth. Protein expression was induced by addition of 0.5 mM 
isopropyl B-D- 1 -thiogalactopyranoside when the optical density at 
600 nm reached a value of 0.6. Cells were allowed to grow for an 
additional 16 h at 28°C. The cells were harvested by centrifuga- 
tion at 8,000 x g for 8 min and the pellet thus obtained was 
resuspended in 40 ml of a solution containing 0.5 M NaCl and 
50 mM Tris-HCl at pH 8.0 (buffer A). Cells were subsequendy 
lysed by the sonication method with an ultrasonic cell-disruptor 
instrument (Tommy) for 15 min (Output 7, Duty 50). A compact 
pellet containing the soluble intracellular components was 
obtained by centrifugation at 40,000 x g for 30 min. The soluble 
fraction was collected and applied onto a Ni-NTA column 
(Novagen) equilibrated with buffer A. Protein was eluted with 
stepwise increase of imidazole (10, 20, 30, 50, 100, 200, and 
300 mM) in buffer A. The fractions containing MCP-1 were 
pooled, and subjected to size exclusion chromatography using a 
HiLoad 26/60 Superdex 75-pg column (GE Healthcare) equili- 
brated with 50 mM Tris, NaCl 500 mM, EDTA 1 mM at pH 7.4. 

Expression and purification of 11K2 scFv 

The DNA encoding the single-chain variable fragment (scFv) of 
1 1 K2 was chemically synthesized by GeneArt and sequence- 
optimized for expression in E. coli (Table SI). The 11K2 scFv 
construct was expressed in vector pUTE [24] displaying a hexa- 
histidine tag at the C-terminus. The DNA sequence was confirmed 
by the dideoxy chain-termination method. 

Cells of E. coli strain BL2 1 (DE3) (Novagen) were transformed 
with the expression vector of 1 1K2 scFv and grown at 28°C in LB 
broth. Protein expression, cell harvesting, and cell lysis were 
performed as described above for MCP-1. A compact pellet 
containing the insoluble intracellular components was obtained by 
centrifugation at 7,500 xg for 30 min. SDS-PAGE analysis and 
western blotting were conducted using the insoluble fraction 
(Figure S2). The soluble fraction was discarded. The insoluble 
fraction was then solubilized with 6 M guanidine-HCl, 0.5 M 
NaCl and 50 mM Tris-HCl overnight at 4°C. After solubilization, 
1 1 K2 scFv was purified in a Ni-NTA column (Novagen) as 
described above for MCP-1, except that the equilibration and 
elution buffers were supplemented with 6 M guanidine-HCl 
(denaturing conditions). 

The purified antibody was refolded by the stepwise dialysis 
method [25]. Briefly, 1 1K2 scFv was diluted to 7.5 uM with 6 M 
Guanidine-HCl in 0.2 M NaCl, 50 mM Tris-HCl, 1 mM EDTA 
(pH 8.0), followed by stepwise dialysis to remove the denaturant. 
In order to increase the refolding efficiency, 0.2 M L-Arg-HCl was 
added to the dialysis solution to minimize protein aggregation 
when the concentration of guanidine-HCl decreased to 1-0.5 M. 
The refolded antibody was further purified on a HiLoad 26/60 
Superdex 75-pg column equilibrated with a solution containing 
0.2 M NaCl, 50 mM Tris-HCl and 1 mM EDTA (pH 8.0). The 
same protocol was employed for the purification of the antibody 
muteins. 
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Figure 1. Crystal structure of 11K2 Fab in complex with its antigen MCP-1. The coordinates were retrieved from the PDB (accession code 
2BDN). (A) Interaction surface of 11K2. The surface in magenta (V H ) and yellow (V L ) correspond to residues of the antibody interacting with the 
antigen. Other residues are shown in light gray. The dotted line indicates the boundary between V H and V L chains. The antigen is depicted in dark 
green. (B) CDR of the antibody. All the residues belonging to the CDR were subjected to virtual systematic mutagenesis. The ribbons in magenta and 
yellow belong to V H and V L chains, respectively. The dark green ribbons depict the antigen. 
doi:1 0.1 371 /journal.pone.0087099.g001 



Binding Assays by SPR 

The interaction between MCP-1 and wild-type 11K2 scFv (or 
muteins) was analyzed by SPR in a Biacore T200 instrument (GE 
Healthcare). Research grade CM5 Biacore sensor chip (GE 
Healthcare) was activated by a short treatment with N-hydro- 
xysuccinimide/N-ethyl-N'-(3-dimethylaminopropyl) carbodiimide 
hydrochloride, followed by immobilization of the antigen MCP-1 
at a surface density of ~220 RU. The activated groups on the 
surface of the sensor were subsequently blocked by injecting 100 ul 
of a solution containing 1 M ethanolamine. The kinetic data of the 
binding of 11K2 scFv to the antigen were obtained by injecting 
increasing concentration of antibody into the sensor chip at a flow 
rate of 30 nl/min. The measurements were carried out in PBS 
containing 0.005% (v/v) Tween-20. Contact time and dissociation 
time were 5 min and 20 min, respectively. Data analysis was 
performed with the BIAevaluation software (GE Healthcare). 
Association (k on ) and dissociation (k 0 fj) rate constants were 
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Figure 2. Mutagenesis in silico. Average energy values (as the sum 
of electrostatic and van der Waals energies) of all possible mutations of 
two different residues (L-Asn31 and L-Ser53) of antibody 11K2 with 
respect to wild-type. Negative values suggest higher affinity between 
the mutated protein and the antigen. The mutations indicated by the 
asterisks were selected for further examination by SPR. 
doi:1 0.1 371 /journal.pone.0087099.g002 



calculated by a global fitting analysis assuming a Langmuir 
binding model and a stoichiometry of (1:1). The dissociation 
constant (Kd) was determined from the ratio of the rate constants 
[26]: 

Ko=k 0 ff /k on 



Calculation of thermodynamic parameters 

Changes in enthalpy {AH 1 ) and entropy (AS") were calculated 
from the slope and intercept, respectively, of the temperature 
dependence of the dissociation constant using the van't Hoff 
approximation [27]: 

In K D = AH" / RT + AS'/R 

where R is the gas constant and T is the absolute temperature. 

The activation energy parameters were obtained from the 
temperature dependence of the association rate constant following 
the Eyring approximation: 

\n(k on /T)= -(AH t /RT) + (AS t /R) + \n(k B /h) 

where k on is the association rate constant, AH f is the activation 
enthalpy, R is the gas constant, T is the absolute temperature, 
is the activation entropy, kg is the Boltzmann's constant, and h is 
the Plank's constant. 

Results 

Computational Selection of Favorable Mutations 

To improve the high-affinity of antibody 1 1K2 for its antigen 
MCP-1 we selected potentially favorable mutations from a virtual 
screening (in silico) consisting of 1,178 single mutations (19 
mutations for each of the 62 residues composing the CDR loops). 
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Figure 3. Binding sensorgrams. (A) Binding of wild-type 1 1 K2 to its antigen MCP-1. (B) Binding of mutein L-N31R to MCP-1. The arrows pointing 
downward indicate injection of running buffer with 11K2 antibody. The arrows pointing upward correspond to the injection of buffer with no 
antibody. The response signal is proportional to the amount of scFv 1 1 K2 binding to a chip decorated with antigen MCP-1 . The straight dotted line at 
the top curve in each panel is drawn to appreciate the slower dissociation rate of the mutein with respect to the wild-type protein. The concentration 
of antibody is given in each panel. 
doi:1 0.1 371 /journal.pone.0087099.g003 



The force field AMBER99 as implemented in the software MOE 
was used to perform 100 energy minimizations of each mutation 
(total was 1 1 7,800 minimizations), and the corresponding values of 
energy were subsequently plotted as histograms (Figure SI). The 
overall shape of the histogram and the median were used to 
estimate the relative efficacy of each mutation with respect to wild- 
type antibody. The relative energies of two representative sets of 
virtual mutations are shown in Figure 2. The residues examined in 
the example are L-Asp31 and L-Ser53, which are located in the 
first and second CDR of the Vl chain, respectively. In a large 
number of virtual mutations, the change of energy is compara- 
tively small (less than ±3 kcal mol -1 ). The greatest differences are 
observed in virtual muteins displaying charged residues (Arg, Lys, 
Glu or Asp). In some cases, the substitution by a charged residue 
increases the attractive energy, whereas in other cases the value of 
energy becomes clearly unfavorable. For example, the calculated 
energy of mutein L-N31R is clearly more advantageous for 
binding the antigen than that of wild-type antibody (A£ L N 51R = 

Table 1 . Kinetic parameters of binding of scFv 1 1 K2 to MCP-1 
at 25°C. 
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— 15.4 kcal mol '), whereas the relative change of energy in 
mutein L-S53R is clearly destabilizing (A£ L 53R = 8.9 kcal 
mol - ). Mutations displaying very favorable changes of energy in 
Figure 2 were selected for experimental validation by the 
technique of SPR (muteins L-D31R, L-D31K, L-S53D, and L- 
S53E). The following muteins were also selected from the whole 
virtual screening and examined by SPR: H-L27R, H-L27K, H- 
N28D, H-N28Q, H-D31E, L-Y30K, L-T56D, L-T56E. We note 
that 1 1 muteins from a total of 1 2 selected mutations involved 
charged residues (92%). 

Evaluation of the affinity by SPR 

The kinetic rate constants of the binding of wild- type scFv 1 1K2 
(and single-muteins) to immobilized antigen MCP-1 were exam- 
ined by SPR (Figure 3, Table 1). The injection of the antibody to a 
surface decorated with antigen produces an increase of the SPR 
signal that is correlated with the association constant rate (k on ) 
(Figure 3A). The dissociation rate constant (k 0 fj) is determined 
from the signal decay after depleting the solutions from antibody. 
The values of k on and k a g determined at 25°C were 14xl0 4 M" 1 
s 1 and 1.0x10 4 s ', respectively, corresponding to a dissoci- 
ation constant (K^,) of 0.8 nM. The small values of k 0 ff indicate 
slow dissociation rates — a clear evidence of tight binding of the 
antibody to the antigen. 

The values of k on , k„ff, and K D were also determined for each 
mutein (Figure 4, Table 1). Significant differences emerge from the 
comparison of their binding affinities. Whereas a majority of 
mutations of the CDR of the Vl chain give rise to a robust 
increase of affinity with respect to the parent antibody (70% of 
mutations), all mutations belonging to the Vh chain destabilize the 
antibody-antigen complex. Overall, the mutation with the most 
favorable effect for the affinity is L-N31R. This mutein binds 4.7- 
fold stronger to the antigen than the wild-type antibody (Figure 
3B;& 0 „ = 13xl0 4 M _1 s" 1 ; k off = 0.22 x 10" 4 s" 1 ; K D = 0.17 
nM). Every mutein exhibiting higher affinity for the antigen than 
that determined for wild-type antibody also displays slower k a ff 
values. In contrast, the destabilizing mutations, without exception, 
accelerate the dissociation of the antibody from the antigen. Thus 
the simple examination of k 0 ff predicts the outcome of the 
mutation in this particular antibody-antigen system. 

Because the mutein L-N31R (displaying the most favorable 
effect on affinity) incorporates the positively charged residue Arg 
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Figure 4. SPR analysis of selected mutations. (A) Location of the residues selected in the virtual screening within the crystal structure of the 
antibody-antigen complex (PDB entry code 2BDN). Mutations belonging to the heavy and light chains are depicted in magenta and yellow, 
respectively. The antigen is shown in dark green. The hot-spot residue PhelOl is also shown (light green). (B) Relative kinetic parameters of the 
binding of the muteins with respect to wild-type protein. Data is given in Table 1. 
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we next examined the effect of the ionic strength in three different 
solutions containing 137, 300 and 500 mM NaCl (Figure 5). No 
major differences are observed in the kinetic rate constants (k on or 
k 0 ff) or the affinity constant (K D ) of wild-type antibody. Similarly, 
the kinetic parameters do not change dramatically in mutein L- 
N31R, although we note that the values of k on decrease slowly but 
progressively from a value of 1 7 x 1 0 M~ 1 s~ 1 in 137 mM NaCl, 
to a value of 13xl0 4 M" 1 s" 1 in 500 mM NaCl (25% decrease). 
Similarly, the affinity also decreases by approximately 25% as 
manifested by an increase of K D from a value of 0.36 nM to a 
value of 0.47 nM. The data indicate a modest role of electrostatic 
interactions, but only during the association (k on ) phase. 

Thermodynamic characterization 

Thermodynamic parameters for the wild-type antibody and for 
the optimized muteins were obtained from the temperature 
dependence of the dissociation constant, as described previously 
(Figure 6, Table 2) [28,29] ._ENREF_28 The van't Hoff enthalpy 
(AH°) and the entropy (—TAS°, calculated at 25°C) corresponding 
to the binding of scFv 11K2 to MCP-1 displayed negative values 
{AH° WT = -7.3kcal mol" 1 , -TAS° WT = -5.0 kcal mol" 1 ) 
indicating favorable contributions from both energetic terms to the 



free energy of binding (AG°wt - — 12.3 kcal mol ). Important- 
ly, the contribution of the enthalpic term increased substantially in 
the muteins. For example, the value of AH° oi L-N31R is 3.5-fold 
more favorable to binding than that of wild- type antibody (AH a L . 
N3iR — —25.6 kcal mol i ,AAH° L . N3JR = —18.3 kcal mol ). 
The change of binding enthalpy of L-N3 1 R is largely (but not 
completely) compensated by unfavorable changes of entropy 
(-TAS° L . N3IR = 12.3 kcal mol" 1 , -TAAS° L . NnR = 17.3 kcal 
mol" 1 , T = 25°C) resultin g in a small advantageous change of 
free energy with respect to the wild- type antibody {AG°l-n31R = 
-13.3 kcal mol" 1 , AAG° L 

-N31R — —1.0 kcal mol ). Similarly, 
the other muteins exhibit favorable changes of enthalpy not 
completely compensated by the entropy term. The thermody- 
namic analysis clearly demonstrates the favorable contribution of 
the enthalpy to the improved affinity, suggesting that the 
mutations generate additional non-covalent interactions between 
the antigen and antibody in agreement with the in silico 
calculations performed above. 

Energetic analysis of the transition state 

The activation energy of each antibody-antigen complex was 
determined from the temperature dependence of k on (Table 3, 




WT L-N31R WT L-N31R WT L-N31R 



Figure 5. Effect of the ionic strength. (A) Association rate constant {k on ), (B) dissociation rate constant (k 0 fft, and (C) dissociation constant (K D ). 
The kinetic parameters of the binding of wild-type 11K2 (or mutein L-N31R) to the antigen MCP-1 were determined in running buffer containing 
three different concentrations of NaCl (137, 300, or 500 mM) at 25°C. 
doi:1 0.1 371 /journal. pone.0087099.g005 
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Figure S3). The activation free energy of wild-type antibody is 
defined by the unfavorable interactions in the transition state (AH* 
> 0), reflecting the negative contribution of the dehydration and/ 
or remodeling of protein-protein interactions during the rate 
determining step (AG^ WTassoc = 10.6 kcal mol" 1 , AH x WT . lss „ c = 
11.4kcal mol" 1 , -TA& WT ,< m oc = ~0.8 kcal mol" 1 ). The 
relative activation free energy of the muteins does not change 
significantly with respect to the wild-type antibody 
(— 0.5<zlzlG : ' : MUT,assoc < 0.3). In contrast, the change of enthalpy 
of the muteins is more advantageous (less unfavorable) than that of 
wild-type antibody (— 6-9<AAH 1l ' MUTilss „ c < — 10.7 kcal mol '), 
suggesting the formation of additional non-covalent interactions 
between the optimized antibody and the antigen during the 
transition state (Figure 7 A). The values of change of enthalpy in 
the transition state are correlated with the values of change of 
enthalpy in equilibrium (Figure 7A). The negative values of 
enthalpy at equilibrium (AAH 0 MU t) and at the transition state 
[A AH* MUT.nssac) demonstrate that the charged residues introduced 
in the optimized antibody improve the enthalpic contribution to 
binding. These observations suggest that the charged residues 
establish electrostatic interactions with the antigen, as depicted in a 
model of the antibody antigen complex (Figure 7B). The novel 
interactions are formed early in the complexation reaction, since 
they play an important role early in the energetic profile of the 
transition state. In the transition state, these non-covalent 
interactions are perfectly counterbalanced by unfavorable changes 

Table 2. Thermodynamic parameters of scFv 11K2. 



of entropy, reflecting the loss of configurational energy at the rate- 
limiting step incurred by the approaching proteins (—TAAS* MU . 
T.assoc = 7.1 ~ 10.7 kcal mol , calculated at 25°C) (Figure S4). 
Although the free energy barrier that the muteins overcome in the 
transition state is nearly identical within experimental error to that 
determined for wild-type antibody (A AG* ~ 0), their energetic 
pathway towards the antibody-antigen complex in equilibrium 
differs from each other. 

Discussion 

The affinity of the therapeutic antibody 1 1 K2 for its antigen, 
the inflammatory cytokine MCP-1, was enhanced up to ~5 fold 
by employing structure-based computational design. Engineered 
versions of the antibody were first designed in silico and 
subsequently verified by in vitro techniques using a recombinant 
scFv construct. The detailed thermodynamic characterization 
revealed the physicochemical principles involved and the operat- 
ing mechanism. 

The selection of suitable candidates of potentially higher affinity 
among a library of 1,178 virtual mutations was carried out by 
computational energy minimizations, and the most promising 
candidates ( 1 2 single-mutations) examined experimentally by SPR. 
Five candidates were mutated at the Vh chain, whereas seven 
candidates were mutated at the Vl chain. In five cases the affinity 
of the optimized antibody increased with respect to the wild-type 
antibody (42% of muteins tested) — a high success rate. 





Protein 


JH° (kcal mol 1 ) 


JJH° (kcal mol 1 ) 


-7i«°(kcal mol 1 ) 


-7^5° (kcal mol 1 ) 


^G°(kcal mol 1 ) 


JJG° (kcal mol 1 ) 


WT 


-7.3 


0 


-5.0 


0 


-12.3 


0 


L-N31R 


-25.6 


-18.3 


12.3 


17.3 


-13.3 


-1.0 


L- S53D 


-13.3 


-6.0 


0.6 


5.6 


-12.7 


-0.4 


L-S53E 


-15.1 


-7.8 


2.4 


7.4 


-12.7 


-0.4 


L-T56D 


-11.7 


-4.4 


-1.1 


3.9 


-12.8 


-0.5 


L-T56E 


-14.5 


-7.2 


1.9 


6.9 


-12.6 


-0.3 



Values of -TAS 0 and AG are given at 25 C. 
doi:1 0.1 371 /journal.pone.0087099.t002 
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Table 3. Activation energy of association of scFv 1 1 K2 to MCP-1 . 



Protein 


z///(kcal mol 1 


Jjrf (kcal mor 1 ) 


-r^lkcal mol 1 ) 


-7"<^tf 7 (kcal mor 1 ) 


^/G 7 (kcal mol 1 ) 


^G^lkcal mol "') 


WT 


11.4 


0 


-0.8 


0 


10.6 


0 


L-N31R 


3.4 


-8.0 


6.8 


7.6 


10.2 


-0.4 


L-S53D 


4.2 


-7.2 


6.4 


7.2 


10.6 


0 


L-S53E 


4.5 


-6.9 


6.4 


7.2 


10.9 


0.3 


L-T56D 


3.8 


-7.6 


6.3 


7.1 


10.1 


-0.5 


L-T56E 


0.7 


-10.7 


9.9 


10.7 


10.6 


0 



Values of -TA5* and AG* are given at 25 C. 
doi:1 0.1 371 /journal.pone.0087099.t003 



Importantly, only mutations made in the Vl chain led to higher 
affinity (five favorable muteins from a total of seven muteins 
examined, i.e. 7 1 % success rate). Such domain bias contrasts with 
other computationally-based optimizations, in which favorable 
muteins are evenly distributed among Vh and Vl chains (reviewed 
in Kuroda et al, 2012) [13]. The deleterious effect of mutations at 
Vh for the binding of antigen reflects the excellent optimization of 
this region in the wild type antibody, both in terms of interaction 
surface (511 A for V H ; but only 257 A 2 for VrJ and the presence 
of a hot-spot residue (Figure 4) [18,30]. In other words, any 
modification of the Vh chain by site-directed mutagenesis disturbs 
the carefully orchestrated interaction surface with the antigen and 
consequently reduces the binding affinity. 

The single-mutations selected from the virtual screening 
incorporate, in 92% of the cases, a charged residue. This 
observation emphasizes the importance of the electrostatic forces 
in the computational optimization. In two previous studies the 
optimized antibodies incorporated multiple mutations in their 
primary sequences [31,32]. For example, the improvement of the 
binding affinity of an anti-epidermal growth factor receptor (10- 



fold) required a triple mutation [31], whereas the optimization of 
the antibody Y0 1 0 1 for the antigen VEGF (9-fold) is achieved by 
incorporating a total of six mutations [32]. Similarly, other 
examples of antibody engineering not employing electrostatic 
optimization also required multiple mutations as described in two 
separate studies (between three and fourteen mutations) [33,34]. 
In contrast, our best design achieves a considerable increase of 
affinity for a matured antibody (5-fold) but requiring only a single- 
mutation, an approach less likely to alter the three-dimensional 
structure of the antibody. We have not examined two or more 
simultaneous mutations of 1 1K2, an approach that could yield an 
optimized antibody displaying even higher affinity than that of the 
single-muteins generated herein. 

The enhanced affinity of the optimized 1 1 K2 antibodies is 
correlated with slower dissociation rate constants (fc 0 ff) rather than 
faster association rate constants (k on ). On the contrary, the 
destabilizing mutations accelerate the dissociation step, indicating 
that k 0 ff is a valid parameter to predict the effect of the mutation 
on the affinity of this particular antibody-antigen complex. A 
previous report suggested that the dissociation step is a first order 




Figure 7. Analysis of the binding enthalpy. (A) Favorable changes of binding enthalpy with respect to wild-type antibody at the transition state 
(empty bar, AAH*) and at equilibrium (filled bar, AAH J. (B) Suggested model of the new interactions formed at the antibody/antigen contact surface 
upon mutation. Residues depicted in yellow and dark green correspond to 1 1 K2 and MCP-1 , respectively. The conformation of the side-chain of the 
mutated residues was modeled from the Dunbrak library of rotamers [44] as implemented in the program Chimera [45] (the most probably rotamer 
was always selected, except in L-Arg31, where the second most probable rotamer was chosen). Because Lys35 of MCP-1 is not interacting with a 
neighboring residue in the crystal structure, the conformation of this residue was also modeled as above. The dotted lines and distances represent 
putative interactions between the mutated residues and the antigen. 
doi:1 0.1 371 /journal.pone.0087099.g007 
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Figure 8. Energetic model of the optimization of an antibody. The diagram represents the enthalpic optimization of an antibody displaying 
higher affinity for its antigen. The enthalpic advantage acquired during the transition state persists in the complex at equilibrium, although it is 
largely counterbalanced by the unfavorable contribution of the change of entropy (entropy/enthalpy compensation). 
doi:1 0.1 371 /journal.pone.0087099.g008 



reaction whose rate is dictated by the strength of short range 
interactions between the proteins (van der Waals forces, hydrogen 
bond, hydrophobic effect, and salt bridges) [35]. Our results 
convincingly demonstrate that introducing a charged residue in 
1 1 K2 increase the affinity for the antigen and slows down k 0 ff. 

The examination of the thermodynamic parameters (Table 2) 
indicates that the optimization of the affinity of the antibody- 
antigen complex is the result of beneficial contributions of 
enthalpic nature originating in the transition state (AAFfi<Q). It 
is important to understand the underlying mechanism in terms of 
the reaction coordinate diagram or "interaction pathway" (Figure 
8) [36-39]. First, antigen and antibody collide with each other 
forming the so-called encounter complex (usee scale) leading to a 
weakly interacting complex in which the hydration of the protein 
surface is not altered significantly (not shown) [36,40] ._EN- 
REF_35 The encounter complex progress towards the transition 
state, although with low efficiency [41]. In the transition state, the 
antigen-antibody partners lose their hydration layers at the 
interaction surfaces, and a reconfiguration of intra- and inter- 
molecular forces takes place. Our results indicate that the 
beneficial effect of the mutations to the enthalpy in the transition 
state (AAH*<0) is also preserved in the equilibrium complex 
(AAH°<0). In other words, the electrostatic interactions generated 
in the optimized muteins are actively contributing to binding early 
in the reaction coordinate, perhaps as early as in the encounter 
complex, which is a state particularly sensitive to long-range 
coulombic forces [36]. The new and favorable non-covalent 
interactions of the optimized muteins are not entirely translated 
into higher affinity of a similar energetic magnitude, since the 
entropy/ enthalpy compensation effect reduces the influence of the 
enthalpy advantage in the free energy [42,43]. Our data indicates 
that ultimately the enthalpy advantage prevails over the entropic 
disadvantage in the final antigen-antibody complex albeit with low 
efficiency since only 5% of the favorable enthalpy is converted in 
useful free energy (Figure S4). 



In summary, we illustrate the benefits of using computational 
methods to design an optimized antibody with higher binding 
affinity for its cognate antigen. In particular, the incorporation of 
charged residues improve the affinity by a combination of 
favorable enthalpic contributions early in the transition state (or 
even earlier in the encounter complex), and slower dissociation 
rates (k 0 fj)- We hope our study will encourage further investiga- 
tions into the fundamental physicochemical basis of protein- 
protein interactions. 

Supporting Information 

Figure SI Representative examples of energy distribu- 
tion histograms corresponding to the interaction be- 
tween 11K2 and MCP-1. The figure shows the histograms of 
the wild type (black dotted), the favorable mutein L-N31R (solid, 
red), and a neutral mutein L-T52S (solid, blue). Muteins displaying 
favorable histograms with respect to wild-type antibody (i.e. shifted 
towards the left in the figure above) were selected for additional 
analysis. 
(TIF) 

Figure S2 SDS-PAGE of insoluble fraction of 11K2 scFv. 

(TIF) 

Figure S3 Decomposition of energy terms along the 
reaction pathway. The data are given in Table 3 of the 
manuscript. In each plot, the three energetic levels correspond to 
the free antibody and antigen, the transition state, and the 
antibody/antigen complex. 
(TIF) 

Figure S4 Entropy/ enthalpy compensation plot. (A) 

Values at equilibrium. (B) Values at the transition state. 
(TIF) 

Table SI DNA and primary sequence of MCP-1 and 
scFv-HK2. 

(PDF) 
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